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ABSTRACT 

The outer regions of disk galaxies show a drop-off in optical and Ho emission, suggesting a star 
formation threshold radius, assumed to owe to a critical surface density below which star formation 
does not take place. Signs of filamentary star formation beyond this threshold radius have been 
observed in individual galaxies in the Ha and recent GALEX surveys have discovered that 30% of 
disk galaxies show UV emission out to 2-3 times the optical radius of the galaxy. We run smooth 
particle hydrodynamics simulations of disk galaxies with constant density extended gas disks to test 
whether over-densities owing to spiral structure in the outer disk can reproduce the observed star 
formation. We indeed find that spiral density waves from the inner disk propagate into the outer gas 
disk and raise local gas regions above the star formation density threshold, yielding features similar to 
those observed. Because the amount of star formation is low, we expect to see little optical emission in 
outer disks, as observed. Our results indicate that XUV disks can be simulated simply by adding an 
extended gas disk with a surface density near the threshold density to an isolated galaxy and evolving 
it with fiducial star formation parameters. 
Subject headings: galaxies: spiral, galaxies: structure, galaxies: evolution, ultraviolet: galaxies 



1. INTRODUCTION 

The most basic description of galactic-scale star for- 
mation is the Schmidt relation, which assumes that the 
star-formati on rate (SFR) i s a power law functi on of the 
gas density (|Schmidtlfl959h . IKennicuttl (|1998l ) showed 
that this relation accurately describes star formation in 
galaxies over a range of 10 5 in gas surface density and 
10 6 in SFR per unit area. However, beyond a critical 
galactocentric radius, SFRs determined from Ho emis- 
sion are truncated relative to the expectation from the 
Kennicutt-Schmidt relation despite t he availability of 
large reserves of low densi t y cold gas (IKennicuttl 119891 : 
iMartin fc Kennicuttl l200lt I Wong fc Blitz) 12002ft . The 



lack of star formation at large radii has primarily been 
attributed to dynamical disk stability, i.e., 
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where k is the epicyclic frequency at a given radius, 
E c (r) is the critical surface density and c s is sound speed 
of the gas, proportional the velocity dispersion of the 
gas ( a = Cs"/" 1 / 2 , where 7 is the ratio of specific 
heats) (|Toomrei fl96l ISpit^eil 11963 lOinrkl Il972^ . The 
effective velocity dispersion of gas in spirals is roughly 
constant across the disk (e .g., Ivan der Kruit fe Shostald 
Il984t lAndersen et al.l 12006ft leaving k as the only radi- 
ally dependent quantity in the critical surface density. 
If velocity curves of spiral galaxies are approximately 
flat, K, and therefore £ c , fall with r, implying that 
at some radius the disk becomes gravitationally stable, 
which may inhibit star formation. In support of this 
notion, several studies have shown correlations between 
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the radius at whic h star formation ceases and the radius 
at which Q ps 2 (jKennicuttJ fl989l: IMartin fc Kennicuttl 
120011 : IWong fc Blitz) I2002L see ISchavd (|2004ft for a com- 
parison of these results). More sophisticated models in- 
cluding the formation of cold molecular gas have also 
been successful at reproducing observed SFR truncations 
(jElmegreen fc Parravanolll994t ISchavdl2004ft . 

Observations indicate the presence of star formation 
beyond the threshold. Faint, isolated HII regions dis- 
tributed in spiral structures ou t to two optical radii have 
been found in face on galaxies ([Ferguson et al.lll998f) and 
in co ld disks in edge-on galaxies (jChristlein fc Zaritskvl 
2008). Simil ar star formatio n has been seen in the 
Milky Way (|Brand et al.ll200Jft an d dwarf galaxies wit h 
gas densities below the threshold (|van Zee et al.l [l997h . 
A system atic survey by the G alaxy Evolution Explorer 
(GALEX: IMartin et all (|2005f )) recently discovered that 
30% of disk galaxies show extended UV emission beyond 
the optical radius, and presumed star f ormation thresh- 
old, of the galaxy (I Thilker et alJ 2007; Gil de Paz et all 
l200alThilker et alJ 2005: Za ritskv fc Christleinl200"7ft . In 
2/3 of these cases, this emission is structured, lying in 
spiral or filamentary patterns (Type I XUV disks). In 
the other 1/3, the emission takes the form of a large 
zone in the outer regions of the galaxy with an enhanced 
UV/o ptical ratio relative to the inner disk (Type II XUV 
disk) (jThilker et aLll2007l hereafter T07). 

To reconcile outer disk star formation with disk stabil- 
ity, it has been proposed that local densities above E c (r) 
in the outer dis k allow star formatio n beyond the trun- 
catio n radii (IKennicuttl 119891: IMartin fc Kennicuttl 
20011: ISchavd 12004 lElmegreen fc Hunter! 120061 : 
Gil de Paz et all 120071) . The low SFR implies a 
low number of O stars producing Ho emission at a 
given time, or none at all, explaining why this star 
formation was undetect ed in the Ha in some galaxies 
(|Gil de Paz et al.l 120071 ) . Most of these disks have 
large amounts of HI in their outer regions and UV 
and Ha complexes are often coincident with local H I 



over-densities ((Ferguson et al.l [l998, T07) supporting 
this picture. Analysis of the star formation profiles of 
these XUV disk galaxies indicates that the azimuthally 
averaged Kennicutt-Schmidt relation could hold beyond 
the dynamical threshold radius (Boiss ier et al.l I2007L 
Hereafter B07). However, studies of star forming clumps 
in the XUV disk of M 83 indicate that traditional star 
formati on laws, includin g the threshold radius, hold 

locally dDong et al.ll2008l). 

Theoretically, lElmegreen fe Hunterl (|2006j ) have shown 
that XUV disks may not be inconsistent with a star for- 
mation density threshold. In their analytical model, gas 
clumping triggered by spiral density waves, radial vari- 
ations in the interstellar medium (ISM) turbulence, and 
gas phase transitions lead to local regions of star forma- 
tion in the extended gas disk. In this paper, we inves- 
tigate XUV disks using smooth particle hydrodynamic 
(SPH) simulations of an isolated disk galaxy with an ex- 
tended gas disk and a simple, commonly employed, star 
formation prescription based on the Kennicutt-Schmidt 
relation that includes a threshold density for star for- 
mation. We demonstrate that, even with these fiducial 
star formation laws, spiral structure propagates from the 
inner disk to the extended HI disk and produces local 
regions of enhanced of gas density which trigger star for- 
mation similar to that seen in observed Type I XUV 
disks. 

2. METHOD 

We used the SPH code Gadget2 (ISpringeHl2005t ) to 
run simulations of galaxies with extended gas disks. We 
incorporate a su b-resolution mu lti-phas e mo del of the 
ISM according to lSpringel fc Hernqui st (2003), which in- 
cludes radiative cooling. The parameters of the multi- 
phase model are £* = 6.5 Gyr, A n — 3000 and T$n = 
3 x 10 8 . The equation of state is determined by these pa- 
rameters. No additional star formation feedback is em- 
ployed. Star formation from the cold phase (total cold 
gas, our simulation does not treat atomic and molecu- 
lar gas separately) follows a Schmidt volume density law 
Psfr oc p^ as with N = 1.5 that is normalized to give 
approximately the same SFR as the Milky Way. A local 
star formation volume density cutoff of 0.004 M Q /pc 3 is 
included to give a surface density thr eshold of ~ 5-10 
M^/pc 2 , consistent with observations (|Kennicuttl 119891 : 
iMartin fc Kennicuttl[200ll ). 

The progenitor galaxy models were constructed follow- 
ing [SjiringeOiriiiT ((200a) ■ The galaxy is roughly analo- 
gous to the Milky Way with a t otal m ass of M 2 oo ~ 10 12 
M© composed of a iHernquistl (|1990f) dark matter halo 
with c = 9, a stellar disk scale length of 3.7 kpc (we use 
h = 0.72 throughout this paper), an exponential gas disk 
of the same scale length. However, we add an approxi- 
mately constant density gas disk that extends to 45 kpc 
and has a surface density of ~ 3.5 M Q /pc 2 . The bary- 
onic mass fraction is 0.04. Forty percent of the baryons 
are in the stellar disk, 20% are in the exponential gas 
disk and 40% are in the constant density gas disk. The 
gas profile of the galaxy at four times in the simulation is 
shown in Figure [T] The galaxy has no bulge and a 0.01% 
M200 mass black hole in its center. The black hole is rep- 
resented by a sink particle that can accrete mass at the 
Bondi rate (|Springel et alj2005l ). but as we are interested 
in the outer disk its presence does not affect our results. 
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Fig. 1. — Gas density profile of the simulated galaxy at the four 
different times shown in Figure[2] Notice how in the first simulation 
snapshot (0.142 Gyr, blue) the extended gas disk has a flat profile, 
while as the simulation goes on and spiral structure develops, the 
outer disk shows peaks that lie above the star formation threshold 
density. From the simulation, the volume density threshold appears 
to correspond to a surface density threshold of approximately ~ 4.5 
Mq/pc 2 which is shown by the dotted line. 



The models have N, 



halo 



1.2 x 10 5 , N s: 



8 x 10 4 



and N gas = 1.6 x 10 5 particles in each component. The 
Gadget softening length of the baryonic particles and 
halo particles are 0.14 kpc and 0.28 kpc, respectively. 

3. RESULTS 

The time evolution of our simulated galaxy is shown in 
Figure [2] The top row shows the stellar surface density, 
the middle row shows gas surface density, and the bot- 
tom row shows the instantaneous SFR surface density, a 
rough indicator of what is observed in the UV. Since the 
instantaneous SFR is shown instead of SFR over the last 
5-200 Myr, this is a lower limit on the UV emission. 

Figure [2] demonstrates that spiral structure from the 
inner disk propagates to the edge of the extended gas 
disk. In the outer disk, peaks of the spiral waves that are 
above the density threshold form stars and the troughs 
that are below the density threshold do not. This gives 
filamentary star formation in the outer disk, as seen in 
Type I XUV disks (T07). Since this is an isolated disk, 
it shows very regular spiral structure, including ring fea- 
tures seen in the second panel of Figure [2j In reality, 
peaks of spiral structure in disks will be influenced by 
interactions with low mass satellites and halo substruc- 
ture, in addition to the dumpiness of the ISM. The 
star formation in this model is not sufficient to create 
a high surface brightness component of stars in the opti- 
cal, even after 4 Gyr. However, very low surface bright- 
ness populations of stars accompany the XUV features. 
Low surface brightness populations of stars are known 
to exist in outer parts of disks accompanying UV so urces 
(|Zaritskv fc Christleinll2007t iGil de Paz et all 12001 . 

4. DISCUSSION 

4.1. Radial profiles 

In the example presented here, we chose a constant 
density radial profile for the extended gas disk. We ran 
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Fig. 2. — Three snapshots from our simulation showing the 
surface density of stellar particles (top), gas particles (middle) and 
SFR (proportional to UV luminosity, bottom). Scales for each are 
shown in the left panel, the color bar for the stellar particle panels 
is converted to projected K-band surface brightness in units of mag 
arcsec -2 , assuming a constant mass to light ratio of 2 Mq/Lq. If 
averaged radially the outer regions have a surface brightness of ~ 
25 mag arcsec . SFR is proportional to UV luminosity, so the 
SFR density range corresponds to a range in UV luminosity range 
of 7 X 10 21 to 7 X 10 29 ergs s" 1 Hz" 1 kpc~ 2 . The panels are 100 
kpc wide. These images illustrate that fiducial star formation laws 
create XUV morphology; notice how the star formation follows 
the peaks in the gas density distribution. The optical emission 
in the outer regions of the galaxy is near or below the surface 
brightness detection limit for most instruments. However, SFRs 
seen in the bottom images are detectable, and are similar to those 
seen in XUV disks (T07, B07). If this simulation were observed, 
spiral structure in the outer disk would be apparent in UV emission 
without accompanying K-band emission, as in Type 1 XUV disks. 

an additional simulation with only an exponential gas 
disk of the same scale length as the stellar disk and found 
that while star formation took place near or slightly be- 
yond the optical radius, the average density of gas at sev- 
eral optical radii is too low to be locally above the star 
formation threshold, so stars did not form there. The 
choice of a constant density extended disk was motivated 
by examples of Type I XUV disks that have fairly con- 
stant densit y extended HI disks (e . g. NGC 4625 NGC 
5055, M83; IBush fc Wilcoti 12001 iThilker et all [20051 : 
lBosmal[l981f) . Our simulation is not intended to repro- 
duce all XUV galaxies, but to illustrate that outer disk 
star formation of the observed morphology is consistent 
with fiducial star formation laws. 

The gas disk we have used for this galaxy model is 
relatively extended (~ 45 kpc radius), to illustrate that 
the size of the disk is not an impediment to spiral struc- 
ture propagating to its edge. The important quantity 
for comparing to observations is the ratio of the sizes of 
the optical and gas disks. Assuming that about 2.5 scale 
lengths of the galaxy would be easily imaged in the op- 
tical, the ratio of the gas to optical disk radii is about 5 
for this simulation, which is realistic when compared to 
XUV disks, for example, M83 (T07). 

4.2. Longevity 

We expect XUV behavior to persist as long as there 
is spiral structure in the outer disk and the peaks in the 



density distribution exceed the threshold. Some XUV 
emission persists throughout the ~ 4 Gyr duration of 
our simulation. The gas density in the outer regions 
is depleted by less than 1% during the simulation and 
therefore XUV behavior will continue for the lifetime of 
the galaxy, as long as this gas is not removed by another 
process (such as ram pressure stripping or tidal interac- 
tions) and spiral perturbations continue. Although XUV 
emission is long-lived, its appearance is time-dependent 
because it follows the spiral perturbations. For example, 
owing to strong spiral behavior at 1.714 Gyr, as shown 
in Figure [5] and Figure [TJ we see rings of UV emission in 
the outer disk. At 3.857 Gyr, there is weak spiral struc- 
ture in the outer disk, leading to only a few spots of UV 
emission. 

In our simulation, the development of spiral structure 
is seeded by nume rical noise arisin g from the discretized 
dark matter halo (jHernq uistj 1 1 993f h However, the simu- 
lation of XUV disks does not require a particular mecha- 
nism for exciting spiral structure. If there is spiral struc- 
ture in inner disks, as observed, our simulation shows 
that it easily propagates to the edge of the outer disk 
and that UV emission tracing that spiral structure can 
occur. 

4.3. The Kennicutt- Schmidt Relation 

B07 average the gas density and SFR density in ellipti- 
cal annuli of 43 nearby galaxies (~14 of which are XUV 
disks) to determine whether the Kennicutt-Schmidt rela- 
tion holds when including UV emission. They find that 
these galaxies fall on the Kennicutt-Schmidt relation and 
extend to surface densities as low as 0.2 M© pc ~ 2 , but 
with considerable scatter. Because our simulation has 
a constant density outer gas density profile at 3.5 M© 
pc ~ 2 , moderated only by spiral density waves, we do 
not probe star formation in such low surface density gas. 
However, in Figure [3] we present a similar analysis of our 
simulation, binning one snapshot in two bin widths, 2.8 
kpc and 8.3 kpc. 

It should be noted that the simulated star formation, 
shown in FigureEJ is normalized slightly below the empir- 
ical Kennicutt-Schmidt law (dotted line) so that the total 
SFR is approximately that of the Milky Way; however, 
this will only affect the vertical normalization and not the 
shape within this plot. In any case, with either the bin- 
ning, our simulation shows evidence of a cutoff, but the 
behavior can be surprisingly varied. When binned finely, 
bins that lie between spiral density waves and there- 
fore below the density threshold, fall off the Kennicutt- 
Schmidt relation, so bins near the threshold density show 
a wide variation in SFR. However, in coarser bins, this 
effect is smoothed as high and low density regions are 
included in the same annulus. Depending on the mor- 
phology of the spiral density perturbation and the bin- 
ning chosen, the variation in SFR near the threshold den- 
sity can be smoothed to the point where these features 
appear to be scatter around a Kennicutt-Schmidt rela- 
tion. We conclude that data in a Kennicutt-Schmidt plot 
can be influenced by such factors as the binning in addi- 
tion to the underlying star formation process and need 
to be interpreted carefully. The Kennicutt-Schmidt plot 
behavior is likely to depend sensitively on the gas dis- 
tribution and more simulations are needed to determine 
whether the B07 results can be reproduced with fiducial 
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Fig. 3. — Kennicutt-Schmidt relation plotted for one snapshot 
of our simulation. The gas density and SFR density are averaged 
in rings at two different bin widths out to 42 kpc. Note that an- 
nuli with a value of lO _6 M0/pc — 2 are upper limits on the star 
formation rate. 

star formation laws. 

5. CONCLUSIONS 

We show that features similar to those seen in Type 
I XUV disks can be produced with fiducial star forma- 
tion prescriptions, including a star formation threshold 
density, if there is gas present at large radii near the 
threshold density. This is a clear illustration of the fact 
that a star formation threshold density does not imply 
the existence of a star formation threshold radius be- 
yond which star form ation ceases. As propo s ed in earlier 
studies and shown bv lElmegreen fc Hunt er (2006), local 
over-densities, caused by spiral density waves in our sim- 
ulation, allow star formation in gas that, when circularly 
averaged, is below or near the classical star formation 
threshold density. 

It is possible that star formation will differ between in- 
ner and outer disks. This could arise owing to different 
properties of the ISM in the inner and outer regions of 



disks, from, for example, a different ionized fraction ow- 
ing to various levels of self-shielding against an ionizing 
UV background. Currently, our simulations do not have 
the resolution to probe the detailed physics of star for- 
mation. Here, we simply show that any simulation whose 
star formation prescription follows a Kennicutt-Schmidt 
relation with a star formation threshold density and has 
a high enough gas density for local over-densities to ex- 
ceed the star formation threshold density in the outer 
disk will exhibit UV emission in the outer disk. 

Understanding the mechanism for generating XUV 
disks has promising implications for probing galaxy 
structure in the outer parts of galaxies. If the distribu- 
tion of H I can be determined based on the observed UV 
emission, this can be used to infer the extended H I prop- 
erties of disks. Outer disk H I and its star formation his- 
tory could also constrain recent gas accretion in galaxies 
and the formation theory of disks. For these reasons, the 
generation of XUV disks warrants further exploration. 
Simulations having different extended gas profiles (e.g. 
exponential, 1/ R) at various gas densities are needed to 
test whether we can create the range of properties seen 
in XUV disks, including Type II XUV dis ks. Upcoming 
extragalactic H I surve ys such as THINGs (|Walter et al.l 
[20051 ) and ALFALFA (|Giovanelli et~aH 120051 ) will pro- 
vide better statistics on the types of H I profiles seen on 
disk galaxies and how they correlate with XUV disk be- 
havior. The role of interactions, particularly flybys, in 
exciting spiral structure and affecting the morphology of 
XUV disks also needs to be explored. Finally, the appli- 
cation of a radiative transfer code to our SPH results to 
co mpare colors of these clumps to samples such as T07 
and lZaritskv fc Christle in (2007) would provide interest- 
ing constraints on star formation such as the initial mass 
function in outer disks. 
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